function y = initial_function1(x)
y = sin(x);
    